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Abstract. The nuclear-cytoplasmic shuttling of NF-kB is characterized by damped 
oscillations of the nuclear concentration with a time period of around 1-2 hours. The 
NF-kB network contains several feedback loops modulating the overall response of 
NF-kB activity. While IkBo; is known to drive and InBs is known to dampen the 
oscillations, the precise role of A20 negative feedback remains to be elucidated. Here we 
propose a model of the NF-kB system focusing on three negative feedback loops (IkBck, 
IkBe and A20) which capture the experimentally observed responses in wild-type 
and knockout cells. We find that A20, like IkBe, efficiently dampens the oscillations 
albeit through a distinct mechanism. In addition, however, we have discovered a new 
functional role of A20 by which it controls the oscillation period of nuclear NF-kB. 
The design based on three nested feedback loops allows independent control of period 
and amplitude decay in the oscillatory response. Based on these results we predict 
that adjusting the expression level of A20, e.g. by siRNA, the period can be changed 
by up to a factor 2. 
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Introduction 

Nuclear Factor-kappa B, NF-kB, is a family of dimeric transcription factors involved 
in a number of important processes such as immune response, cellular growth and 
apoptosis pQ. NF-kB regulates the expression of more than a hundred genes and is 
implicated in a large number of diseases, including cancer, heart diseases and asthma pQ. 
Nuclear translocation of NF-kB, necessary for its activity, is triggered by a wide variety 
of stress signals: endotoxin LPS, cytokines IL-1 and the tumor necrosis factor (TNF). 
Fluorescence imaging of the TNF-triggered NF-kB activity in single mammalian cells 
shows distinct "spiky" but asynchronous oscillations in the level of nuclear NF-kB [2]; 
populations of mouse fibroblast cells continuously exposed to TNF exhibit damped and 
smooth - probably due to population averaging - oscillations in the nuclear NF-kB 
concentration. The production of damped oscillations with a time period of around 1.5 
hours thus seems to be a robust characteristic of the NF-kB system. NF-kB is regulated 
by several negative feedback loops: two acting through the inhibitor proteins I/tBa, e 
which bind and sequester it in the cytoplasm. Another feedback regulates concentrations 
of nuclear NF-kB through A20, see Fig. [TJ Addition of TNF activates the IkB kinase 
(IKK) which in turn causes the phosphorylation, and subsequent degradation, of the 
IkB inhibitor proteins, thus releasing NF-kB. Free NF-kB translocates to the nucleus 
inducing transcription of the inhibitor proteins, \nBe and IkB«, and A20. In turn the 
I/tB proteins inhibit the NF-kB transcription factor by actively exporting it out of the 
nucleus. A20 on the other hand acts upstream by inactivating IKK, see Fig. [TA. There 
is another inhibitor protein, IkB/3, but it is only slightly induced by NF-kB compared 
to the two other inhibitor proteins [3] resulting in a weak negative feedback on NF-kB. 
We will thus omit it in our model, see Fig. [T]B. 

Each IkB protein forms a negative feedback loop as they are all transcriptionally 
activated by NF-kB. The loops are not identical. Knockout mutant studies in bulk 
indicate that in the absence of InBe the nuclear NF-kB oscillations are enhanced whereas 
there are no oscillations in the absence of InBa. It has been suggested that InBe 
dampens oscillations generated by luBa, which is consistent with the differences in 
their half-lives - InBe is at least twice as stable as InBa [I]. InBa is activated almost 
instantly while InBe activation occurs 37 min after NF-kB enters the nucleus. This 
difference in half life and time of activation allows us to define InBa as a "fast" and 
IkBe as a "slow" negative feedback, see Fig. CD 

A20 feedback acts upstream of NF-kB and IkB. It is an important regulator of 
late IKK activity and was shown experimentally to be required for the drop in NF-kB 
activity separating early and late phase response to TNF when measured in bulk pEj. 
Cells deficient in A20 show persistent IKK activity and develop severe inflammation 
and cachexia [5]. 

The physiological importance of NF-kB transcription factor and its intriguing 
dynamical behavior made it a center of attention for decades both from an experimental 
and theoretical point of view [6l [7]. The first computational model of the NF-kB 
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Figure 1. Schematic drawing of the core of NF-kB regulatory network. A) Details 
of the NF-kB regulation. TNF activates the IkB kinase (IKK) which in turn causes 
the phosphorylation, and subsequent degradation, of the IkB inhibitor proteins, thus 
releasing NF-kB. Free NF-kB translocates to the nucleus inducing transcription of the 
inhibitor proteins, IftBe and IkBc*, and A20. The IkB proteins inhibit the NF-kB 
transcription factor by actively exporting it out of the nucleus. A20 acts upstream by 
inactivating IKK. The IkB/3 is not shown. B) Nested negative feedback perspective on 
the NF-kB regulation. NF-kB is regulated by two parallel negative feedbacks acting 
through the inhibitor proteins IkBo, e and A20 negative feedback, acting upstream by 
inactivating IKK. The variation in the greyscale of the three feedback loops indicates 
the difference in timescales: dark grey stands for fast and light grey for slowest. 



pathway was proposed in Hoffmann et. al [8] and used to understand the dynamical 
responses of the NF-kB wild-type and IkB knockout, e.g. oscillations and their absence 
in knockouts. This model has later been modified and used by Nelson et al [2] to 
analyze oscillations in single cells. Krishna et al 2006 |9] showed that the model can be 
significantly reduced while still capturing the essential dynamical features, in particular 
showing spiky oscillations in single cells. 

Both modeling and experimental results suggest that A20 is important for lowering 
the level of nuclear NF-fcB after TNF stimulus [3 H]. Other studies have focused on 
details of where and how A20 acts in the pathway [10J. To address the discrepancy 
between bulk and single cell data Ref. [TTJ [10] introduced stochasticity and showed that 
averaging single cell stochastic dynamics leads to a smooth damped response in bulk. 

In the current view of the system, InQa feedback drives the oscillations and the 
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Figure 2. Model simulations of the four states of the cell obtained by numerical 
integration of the Ordinary Differential Equations describing the model in Methods 
section. A: wild-type, B: IkBe~/~, C: A20 _/ ~, D: A20~/~ /IkBe~/~ . In each panel 
blue line is the concentration of active IKK, red is the normalized concentration of the 
nuclear NF-kB and dashed green line is the concentration of total inhibitor proteins, It 
= InBa + IkBe. Note the differences in NF-kB oscillation period and decrease in the 
amplitude in absence of A20 and IkBe~/~ and increased levels of IKK in late stage 
for A20 knockout cells. 




2 4 6 8 10 12 

time, hours 



Figure 3. A20 changes the period of nuclear NF-kB oscillations. WT is in black, 
reproduced from Figure 2A, and A20 knockout is in red, reproduced from Figure 2C. 
Note the difference in steady states after oscillations are damped. 



\nS>e feedback dampens the oscillations. From the bulk experiments A20 is known to 
lower the level of nuclear NF-kB after stimulus. It is not known how the temporal profile 
of NF-kB in single cells is affected by A20 [I]. We would like to investigate whether 
and how A20 is modifying NF-kB oscillatory behavior in single cells. To address this 
question we will extend the NF-kB model by Krishna et al. [9] to include IttBe and A20 
negative feedbacks. 
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Model 

Figure Q] shows schematic representation of the model. It contains three negative 
feedback loops centered around NF-fcB: IfcBa, \kBe and A20. 

The dynamical variable of most importance is N n , the nuclear NF-kB concentration. 
The first term in the equation for N n is the rate of increase in nuclear NF-kB 
concentration due to import of free NF-kB from the cytoplasm. This rate is lower 
for higher levels of the IkB proteins. The other two negative terms model the decrease 
of the nuclear concentration due to sequestration by the IkBs and subsequent export into 
the cytoplasm. Over the timescales we are interested in there is no significant production 
or degradation of NF-kB. The mRNA levels of \nBa and \kBe are regulated through a 
sigmoidal function of NF-kB, given by N 2+ K 2 ■ Here we assumed that there is a weak 
cooperativity in NF-kB activating transcription of luBa/e with Hill coefficient two. We 
also include a small basal level of transcription in luBa, which has little influence on 
NF-kB dynamics but has an important role in reproducing the InBa and IkBe mRNA 
fold induction experimentally measured in [3]. 

At the protein level, the rate of protein increase is linearly proportional to the 
respective mRNA. The rate of decrease in luBa is controlled by IKK-independent 
degradation, I a /r a and IKK-dependent degradation. The rate of IuBa decay is 
proportional to both IKK and the concentration of complexes formed between InBa and 
cytoplasmic NF-kB, N c , oc IKK[I nBa : N c ]. Assuming that reaction rates for complex 
formation are much faster then the nuclear import/export and luBa degradation, the 
concentration of [InBa:N c ] can be derived to be (1 — N n )/(1 + IkBe + IuBa) (see 
Supplementary Materials for more details). The changes in the IkBe protein levels 
are governed by the same terms as for luBa, but with different rate constants, which 
overall make the negative feedback through IkBe slower and weaker (see Supplementary 
Material for details). Apart from this, the equations for IkBe are the same as for InBa. 

The IkB Kinase (IKK) is the driving force of the system as its activation leads to 
the degradation of the inhibitor proteins, the IkBs, and thereby the release of NF-kB. 
IKK is activated upon stimulation of the membrane receptor but the detailed mechanism 
for this activation remains to be clarified. We have chosen to use the mechanism earlier 
proposed in Ref. pU] to model the IKK activation: a three step process where IKK 
is converted from its neutral state to being active by the triggering signal, TNF, see 
Fig. [T]A The active IkB kinase can turn itself off and go back into neutral state before 
being activated again by the TNF signal. The IkB kinase is shut down by A20 which 
inhibits the transformation from inactive to neutral IKK thereby leaving IKK in an 
inactive state, see Fig. [TA. 

The first 30 minutes of IKK adaption-like temporal profile in response to TNF 
stimulation - which appears to be independent of A20 regulation - is modeled such that 
IKK peaks after about 15 minutes and then goes to a new steady state after 30 minutes 
of TNF induction. The new steady state is however determined by A20: it is high in 
the absence of A20 and decreases with increasing concentrations of A20. This is a slow 
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Figure 4. Analysis of the lnBe~l~ and A20 feedback loops. The strength of each loop 
is varied by varying the IkBs transcription and the A20 coupling to IKK. The white 
circle represents the reference state of wild type (Figure 2A). The colorcoded are: A) 
The average period taken between the 2nd and 3rd peak. B) The amplitude of NF-kB 
averaged over 2-nd and 3rd peak with a lower boundary of 1%. C) The dampening 
effect, tells how fast are oscillations damped, measured by difference in amplitude 
between the 2nd and 4th peak. Bifurcation line is shown in black and separates 
regions of sustained oscillations from regions of damped oscillations. This bifurcation 
line has been found by counting the number of peaks in the nuclear concentration of 
NF-kB within the simulated timespan (see Methods section for details). Both InBe 
transcription and A20 coupling to IKK is able to change the amplitude (B) and dampen 
the oscillations (C) but only the A20 coupling to IKK is able to alter the period of the 
oscillations (A). 




feedback as IKK must first activate NF-kB leading to the production of A20 that in 
turn shuts down the pathway. 

We have taken most rates and timescales from existing literature wherever possible 
and manually adjusted so that the model reproduces the following experimental 
observations: 

1. Wild-type cells show damped oscillations in nuclear NF-kB with a time period of 
60-100 min. 

2. Mutants with luBa alone show enhanced oscillations with same period of 60-100 
min, while those with InBe alone do not show oscillations. 

3. The fold induction in mRNA of InBe reaches twice the level of fold induction of 
the InBa mRNA 

Results and discussion 

Model validation 

The basic response of our model, with the default parameters, to a continuous presence 
of TNF is damped oscillations of nuclear NF-kB. The original wild type response as well 
as luBe and InBa knockout matches the experimental observations (see Fig. [2A, B), 
Supplementary Fig. SI and interactive applet [12J. The time period of the oscillations is 
about 100 minutes, which also matches experiments. Thus, the basic response (criterion 
1 and 2) is correctly reproduced by the model, Fig. [2j We also checked the model against 
criteria 3 and other known knockouts, for further details see Supplementary material. 
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Figure 5. Set of characteristic responses in nuclear NF-kB obtained at different InBe 
transcription rates and A20 coupling strengths to IKK. The colorcoded is the average 
period taken between the 2nd and 3rd peak (a section of A in Fig. 2]). Several parameter 
sets have been selected both within and outside of the bifurcation border (the black 
line). Both IkBe transcription and A20 coupling to IKK is able to change the amplitude 
and dampening of the nuclear transcription factor activity but only the A20 coupling 
to IKK is able to alter the period of the oscillations. 



It has been suggested that A20 lowers nuclear NF-kB levels of the late phase of the 
response. This observation stems from bulk experiments [131 SJ- However, the role of 
A20 has not been investigated in single cells where NF-kB has oscillatory response [10J. 
Using the proposed model we aimed to investigate the role of A20 as a modifier of the 
NF-kB oscillatory behavior. We find that A20 is able to adjust the period of nuclear 
NF-kB oscillations in the range from 1 to 3 hours. Below we describe the details of our 
finding. 

We have modeled the NF-kB response in A20 knockout cells and were surprised to 
find that not only is the resting level of the late phase of the NF-kB response increased 
(as shown in Fig. S2) but also the oscillations become more pronounced when compared 
to wild-type cells, see Fig. |2] A,C and Fig. [31 The amplitude of the NF-kB response 
is decreased in the presence of A20 whereas the period is increased. This effect of 
A20 dampening the oscillations and increasing the period is even more visible in IkB^ 
knockout cells where the oscillations are sustained, compare Fig. [2J3 and D. 

A20 only has an effect on the temporal profile of IKK in its late phase: the A20 
mRNA level peaks at 30 min and the protein shows its effect after 45-60 min [T3l H]. 
This feature is captured by the model where the level of IKK in the late phase is pushed 
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down in the presence of A20 and generates low frequency NF-kB oscillations. In the 
absence of A20 the late phase of IKK stays at a high level and generates high frequency 
oscillations, compare Fig. [2A and C and Fig. [3j 



IkBe does not change the period of NF-kB oscillations 

Is A20 feedback loop dispensable and can IkBe have similar effect on the oscillations 
period? There are two possible mechanisms for IkBe to dampen luBa oscillations: 
a) destructive interference, where InBe oscillates on the same time scale as IuBa but 
with the shifted phase compared to IkBoj, such phase shift or delay has been observed 
experimentally [3| [TO] and b) where IkBe varies on a slower time scale and the interference 
between fast changing IkBo and slow InBe results in damped oscillations. In order for 
a) to be true InBe and InBa must exhibit comparable strength and frequency and thus 
IuBa knockout should oscillate in single cells, however at present there is no such data 
available and bulk experiments show no oscillations in the InBa knockout. On the other 
hand it is known that InBe has a slower degradation through IKK dependent mechanism 
[H] which makes InBe change on a slower time-scale. Furthermore the b) scenario can 
easily reproduce bulk data for InBa knockout and therefore we model InBe through 
scenario b). With this constraint InBe does not influence the frequency in our model 
(compare [2JA. and B). 

Function of A 20 and InBe is robust to variation in parameter values 

Our observation that A20 decreases the frequency of NF-kB oscillations whereas IkBe is 
not capable of the same effect, is based on the model with specific set of parameters that 
were chosen to fit experimental results. To check if our observation is not merely a result 
of a specific parameter combination we have investigated how the period, amplitude as 
well as the oscillation dampening change with varying some of the model parameters. In 
Figure 0A we show how the period changes as the strengths of A20 and InBe feedback 
loops vary 100 fold above and below the values we used in Figure |2j The period is almost 
independent of the strength of InBe whereas it gradually increases as the strength A20 
feedback loop increases, thus supporting our conclusion that A20 is a key regulator 
of the frequency of NF-kB oscillations. It is important to note that although it is 
possible to decrease the frequency by increasing the strength of the IkBe loop when 
IkBe transcription rate is larger than 50 (see top left corner in Figure HA), the resulting 
oscillations are strongly damped as illustrated in Figure HP. 



Reasons for the differential effect of A20 and IkBe on the period 

We can understand why A20 and not IkBe feedback can effectively adjust the frequency 
by looking at the IKK dependent degradation of IkBo and \kBe terms: 

ikBcx degradation: ~ A- r— ; r— r 

1 J l + [/ a ] + [/ e ] 
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and 

IkBe degradation : ~ [K] - ^ 



1 + [h] + [J e ] 

IKK acts similarly on both IkBe and luBa feedbacks: Higher levels of IKK 
correspond to faster degradation of luBe/a and thus shorter period. In case of IkBe, 
however, - because of competitive binding to cytoplasmic NF-kB - higher levels of 1k.Be 
will lead to faster degradation of IkBe but slower degradation of IkB«, see equations 
above. Thus IkBe can only increase the period when oscillations are driven by I^Ba, 
whereas IKK can work in both directions. There are, however, limits to how much IKK 
can increase the period. This limit is set by the fact that as A20 increases, and IKK 
decreases, the system passes through a Hopf bifurcation (shown by black line in Figure 
HJ) where it goes from sustained to damped oscillations so that the amplitude and period 
of oscillations are strongly correlated: the longer the period the smaller the amplitude, 
compare Figure 0A and B. 



Conclusions 

Growing evidence indicates that temporal control of NF-kB and the downstream genes 
is of crucial importance for cell functioning: constitutively active NF-kB is a cause of 
many human tumors. Active NF-kB turns on the expression of genes that keep the 
cell proliferating and protect the cell from conditions that would otherwise cause it to 
die via apoptosis. At the same time defects inactivating NF-kB result in increased 
susceptibility to apoptosis leading to increased cell death. It appears that the original 
solution to this dilemma is through transient activation of NF-kB [8] which on a single- 
cell level presents as damped oscillations. Such temporal control can allow for selective 
gene activation [T3l [TO] . Given that the NF-kB temporal response is of a high importance 
and is primarily regulated by three negative feedback loops we investigated the role of 
each negative feedback in shaping the response with the main focus on A20 negative 
feedback. 

We find that the design of having two luB feedback loops on the same level as 
NF-kB and the upstream feedback from A20 allows for a wide variety in possible NF- 
kB outputs. The fast feedback from IuBa generates the oscillatory behavior of the 
transcription factor which is damped by the delayed and out of phase activity of IkBe. 
This general response of NF-kB upon TNF stimulation is altered by the upstream 
feedback from A20 which is slower than the two other feedbacks and acts on the IkB 
kinase and not on the transcription factor itself. By changing the IKK profile, A20 is 
able to modulate the frequency of the transcription factor in a way not possible from 
the two IkB feedbacks. Additionally, A20 lowers the end level of the nuclear NF-kB 
meaning that when the system comes to rest more NF-kB is removed from the nucleus, 
in the presence of A20, compare Fig. [2] A and C and Fig. [3] [4J. (This is because IKK 
is less active and thus not degrading the IkBs as fast, leaving the inhibitor proteins to 
enter the nucleus and export the transcription factor). 
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If the biphasic response - as seen in bulk experiments - is a result of the population 
average of single cells with oscillating NF-kB, then, in bulk experiments, A20 will exhibit 
its effect by affecting the timing of the second phase onset. Thus a specific prediction 
would be that - as the timing between first two peaks is shorter, see Supplementary 
Figure S3, - the second phase should start earlier in A20 knockout cells. 

We found that similarly to IkBe, A20 is able to dampen oscillations albeit through 
a different mechanism. Whereas IkBe dampens oscillations through competition with 
\nBa for binding to NF-kB and acts about an hour after TNF induction, A20 dampens 
oscillations by inhibiting active IKK during late phase of the response, after about 
2-3 hours. A distinct novel finding of our investigation suggests that not only does 
A20 dampen the oscillations, it can also adjusts the period of the oscillations. Thus, 
A20 together with I/tBe allows independent tuning of both dampening - through 
\kBe negative feedback - and the frequency - through A20 negative feedback. This 
combination of nested feedback loops covers a wider variety of temporal responses where 
one can access both sustained oscillations and damped oscillations with low or high 
frequency. 

Our findings lead to a clear prediction that in single cells decreasing the coupling 
between A20 and IKK should lead to higher frequency oscillations in NF-kB. This can 
be experimentally tested by knocking down A20 with siRNA. Our model also predicts 
that the frequency of oscillations can be further increased by increasing the TNF dose in 
A20 knockout cells, where the IKK level in the late phase of the response is proportional 
to TNF dose. An interesting future direction would be to examine how the diversity 
of NF-kB oscillating temporal profiles created by nested feedback loops can allow for 
selective gene activation. 
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Methods 

Model Description 

dN n _ (N t -N n ) N n 

dt A (K I + V l a + V l e ) * [Ia + l£) (8 + N n ) [A) 

I N 2 

~dT~ V a {K 2 D + Nl) lmama { ' 

dl a v (N t -N n )I a 

dt {Kj + rjl a + rjl £ ) 

dime A„ 
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dA , , 



(.5) 



(.6) 




(•7) 



— = T(K t -K- Ki) - fiK 2 




(.8) 



(.9) 



(.10) 



The equations are all rescaled, see Supplementary Materials for details of the rescal- 
ing and the parameter values. 



Approximating Bifurcation Line 

To find the border for the bifurcation we counted the number of peaks in the nuclear 
concentration of NF-kB. We counted the oscillations as sustained if the number of peaks 
was 100 or above. This approximation slightly overestimates the size of the region with 
sustained oscillations and is not affecting any results and conclusions. 
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Table 1. Model variables and parameters 
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Figure 1. Model simulations of the four states of the cell : A: wild-type, B: IkBe 
C: A2Q-/-, D: A2Q~>~ /IkBe'I- 
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Figure 2. Model simulating the wild- type cell. Left column: Recordings of the 
durations of the I s *, 2 nd and 3 rd period measures as the time between first and second 
peak, second and third and trird and fourth peak. Right column: Amplitude of I s *, 
2 nd and 3 rd peak. 



Initial Model for IkBs regulation 

We use the following abbreviations: N n &N, free nuclear and cytoplasmic NF-kB; 
I m , IkB mRNA; free nuclear and cytoplasmic IkB; (NI) n &c(NI), nuclear and 

cytoplasmic NF-kB-IkB complex; IKK, IkB kinase. 

The equations for IkBe are the same as for InBa so we just use variable I in our 
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derivations. 

The seven- variable model is defined by the equations 

km n N - k fn N n I n + kbn(NI) n , 



dt 



dIm _ k N 2 _ j 

" J t ly n lm 1 mi 



dt 

^- = k u Im - kfNI + k b (NI) - k Iin I + k Iout I n , 
dt 

dN 

— = -k f NI + (k b + a)(NI) - k Nm N, 
= k f NI - (k b + a)(NI) + k NIout (NI) n , 

—jj- = kn n I — ki ou tI n — kf n N n I n + k bn (NI) n , 

= kfnN n I n - (k bn + k NIout )(NI) n . 
Reduced Model for IkBs regulation 

First, taking note of the fact that kf and kf n are large, we assume that all complexes 
are in equilibrium, i.e. 

kfNI « (k b + a)(NI), 

kf n N n I n « (k bn + k NIout )(NI) n . 

Simulations show that these are good approximations. In terms of J*°* = I n + (NI) n 
and N l c ot = N + (NI) = N tot — N n , which are slowly varying, we can rewrite the above 
equations as follows: 

(NI) = (N tot - N n ) 



Ki + V 



N = (N tot -N n ) 1 ^ n , 
(NI) n = I M - V " 



K N + N n ' 

j jtot K N 

n n K N + NJ 

where Kj = (k b + a)/kf = 0.035 [iM and K N = (k bn + kNiout)/kf n = 0.029 /iM, using 
the parameter values above. 

Using these expressions, the equations of the seven-variable model reduce to the 
following four (Fig. 9): 

dN n _ (Ntot - Nn) , I™N n 



dt Ki + I K N + N n 



Theoretical analyses predict A20 regulates period of NF-kB oscillation 17 

dim 



k t Nl - 7 m i m , 



dt 



dl (Ntot ~ N n )I I™ 
- ktUm - a — kj in l + k Iout K N - 



rft A/ + J A/v + iV n 

n — h T h K n h n n 

^Iin 1 r^Iout-^NT^ ; TT - &NIout~ 



at K N + JM n K N + JM n 

jtot 

* N K N +N n 



jtot 

First, we note that the terms —ku n I and ki out K^ K " N in the dl /at equation 



are much smaller than ~ ^ Nt K I +f >I an d can be neglected as long as IKK is nonzero. 
Second, simulations reveal that the term kwiout k" n +n ' * n ^ ne dl^/dt equation, also 
shows sharp spikes as a function of time which coincide with the spikes of N n . The 
value of this term is substantial only when N n 3> K N , i.e., during the spikes of N n , and 
at those times J*°* dips to its minimum. We therefore make the approximation that 
can be replaced by its minimum value, I^ min , which satisfies the equation 

jtot jtot J\J 

1 j 1 ry- n,min . i n,min " 

Klin 1 — Kiout -ft- TV 77 ; — r I^NIouf 



K N + N n K N + N n 

In the regime where N n ^> K n this gives 

jtot ^ _klm_ 
n,min ~ 7 

K NIout 

Using this we can reduce to a three- variable model 

ciiV n (JVg - N n ) _ IN n 



dt ^ m 1 Kj + I " lm 5 + N n ' 

dim , KT 2 j 

dl , r (A^ - N n )I 

Tt = ktlIm ~ a k i + i ■ 

Parameters and variables re-scaling 

We start with the following system of equations 

dN n (N, - N n ) N n 

~dt~ = h ~ Kl w, + t. + l) - B{1 - + ^wrm <■"> 

I N 2 

1 ma _ , iv n j I 1 r>\ 

dt ~ P+ a {Kl + Nl) lma ma { ' 

dl a ^ (N t - N n )I a I a 

-7T = k Ua I ma - a a K — — (.13) 

at A/ + i a + i £ 7 a 

_ . N n I m£ , s 

dt " re (^ + iV2) 7me l ' } 

dl e T ^ (N t - N n )I e I e 

— = ktlelma ~ Ci e K———— (.15) 

at Kj + I a + I £ 7 e 
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= t A Nl/K DA - (.16) 

dt ~lAm 

dA , A A , -, ^ 

— = fctuAn (-17) 

a* 7a 

a — = T{K t -K-K l )-iiK' 1 (.18) 

%=^-0Jri ( - 19) 

(.20) 

where fcj m is the import of the NfkB into the nucleus, Kj is dissociation constant 
for IkBs binding to Nfkb, ... 
Re-define 

1 Re-define N n by scaling with Ko and N t 

N' n = N n /K D , N' t = N t /K D , 5' = S/K D , a' a = a a /K D , a' e = a £ /K D 
and furthermore, 

K = K/K = N n/N t , 6" = 5'/N' t = 5/N t 

After this transformations and redefining N" = N n , 5" = 5, etc.the equations will 
be rescaled to : 

dN n T ^ (1 - N n ) N n 

— = hmK! - B(I a + I £ )- 



dt im <(K I + I a + I £ ) v a £J (5 + N n ) 
I N 2 

1 ma , Jv n _ j 

j. P ~r laTj ' , r9 x Ima^ma 

dt (1 + N*) 

dla_ k r K (1 ~ N n )I a I* 

dt K j + I a + 4 7„ 

the equations for I £ and I m£ scaled similarly and the rest remains as above. 

2 Re-define I a and I £ by scaling with Ki. 

U = ^^B' = BK T , k' tla/e = k' tla/ jKj, a' a/£ = a a/£ /Kj 

3 Re-define I am and I em by scaling with k' tla . 

Re-defining I' a/em = k' tla/e I a/em = k tta i e / Kil a/em , p' = pk Ua/e / K u t' aje 
ta/ektia/e/ Ki results in following equations for 

dNn -k m ^^-B { i a + i £ y 



dt ™(1 + J a + J e ) V a ^ {6 + Nn) 



I N 

1 ma 1 ' • 



2 



dt a a+ND 

dla _ (l-N n )I a _ I« 

dt 1 + I a + I £ 7 Q 
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4 Re-define A m by scaling with k t \A- 

A' m = k t iAA m , =>■ t' A = kuatAl Kda, °~ = °~/ku a 2 thus equations for A20 become 

dA A 

= t A iy n 



dt n 'jAm 

dA A 

at 7a 
5 Re-define K and Ki by scaling with K t . 

K = K/K t , K,i = Ki/ K u jj, = K t fi thus equations for IKK become 

^=T(1-K-K t )-^ 



dt ^ r a A 2 + 1 



Thus the final system of equations is: 



dt "™ (1 + I a + I e ) ~^ a '^(5 + N n ) 
I N 2 



dt * °(1 + JV*) 

dIa_ T _ T ^ (l-N n )I a I a 

j — J-ma Ct a J^ i , t , t 

dt 1 + I a + I e 7 a 

dt le (l + N 2 ) lmeme 
dh T ^ (l-N n )I £ h 

dt ~ lme aeK l + I a + I e l£ 

dAm - t A N 2 n - Am 



dt n ^Arn 

dA A 

—77 — A rn 

dt 7a 

dK , , 

— = T{l-K-K i )-nK 2 

dt P 1 aA 2 + l 



.21) 
.22) 
.23) 
.24) 
.25) 
.26) 
.27) 

.28) 

.29) 
.30) 



and the scaled variables and parameters are: 

AT _ Nn. J — J — k tla/e. J A - U, . A ■ K — —■ K ■ — ^ 

n Nt ' a / £ Ki ' ma /£ Kj ma/si ^-m <HlA- r *mi J1 j( t i sy i 

5 = £ = 0.0414; B = BKj = 0.014; fc« a/c = %f ;; p = = 58.4; 

*a/e = t a/e ^,t a = 7300, t e = 27.4; < /£ = = 219, a e = 6.57; 

t A = ktu-S*- = 18; a = = 77xl0" 5 ; (3 = 5; u = iQit = 100; TNF changes from 

^DA tlA 

0.001 to 2.5; = 6.67 
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Variable 


Description 




nuclear NF-kB normalized to total NF-kB 


±a / e K l 


free IkBs scaled with dissociation constant Kj of I/tBs binding to NF-kB 


J _ ^tla/e J 
^ma/e J-ma/e 


re -defined value of IkB mRNA, ku a / e is the translation rate of IkBoc/e 


A m — ktlAA rn 


re-defined A20 mRNA, kuA is the A20 translation rate 


K = ^ 

Kt 


active IKK normalized to the total IKK, K t 


K = & 

ts.t. 


inactive IKK normalized to the total IKK, K t 


Scaled Parameter 


Description 


6 = ± (fiM~ r ) 


concentration at which half of the IkBa/e is bound in complex with 
NF-kB, normalized to total NF-kB 


B = BKj 


proportionality factor of the export of nuclear NF-/tB, 

scaled with the respective translation rates and dissoc. constant of 

NF-kB binding to IkBs, Ki 


A = AKj 


proportionality factor of the import of NF-kB, 

scaled with the respective translation rates and dissoc. constant of 

NF-kB binding to IkBs, if/ 




constituitve, NF-kB dependent transcription rate 


of IkBa mRNA, scaled with the respective translation rates and Ki 


4- ..— + ktla/e 
ta/e ^a/e 


NF-kB dependent transcription rates of IkBs mRNA scaled with Kj 


tA - hlA K ^ A 


A20 transciption rate scaled with A20 translation rate, kuA, 


and diss, constant of NF-kB binding to DNA at the operator site 
controlling A20 promoter. 


llrnal I Ima 


half-life of \kBci mRNA scaled with 7/ ma 


lime/ lima 


half-life of I/tBe mRNA scaled with 7/ ma 


1 1 a/ el Ima 


half-life of the I/tB's scaled with 7/ ma 


lA20m1lma 


half-life of the A20 mRNA scaled with 7/ ma 


lA20llma 


half-life of the A20 scaled with 7/ ma 


aa / £ K D Kj 


rate constant for IKK dependent degradation 
scaled with dissoc. constant of NF-kB binding to 
operator site at the IkB's promoter, K D and Kj 


\x = K t fi 


rate of IKK self-inactivation scaled with total IKK, K t 


K tlA 


strength of A20 negative feedback scaled 




with the square of A20 translation rate, k\ lA 



